Neutron scattering evidence on the nature of the boson peak 
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A close inspection of neutron spectra of glass formers in the frequency region above the boson peak 
does not show the constant eigenvalue density expected for a random dynamical matrix, but a slow 
decrease towards higher eigenvalues. In 1,4-polybutadiene and selenium, the slope becomes steeper 
with increasing temperature in the undercooled liquid, suggesting an influence of the vibrational 
entropy. One can describe the behaviour quantitatively in terms of a balance between level repulsion 
(from the randomness of the dynamical matrix) and vibrational entropy. On the basis of this 
evidence, the boson peak must be a crossover from such a balance at higher frequencies to a mixture 
of sound waves and additional excitations at low frequency. 

PACS numbers: 63.50.+X, 64.70.Pf 
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There is as yet no generally accepted explanation of the 
bosonpeak in the neutron or Raman spectrum of glasses 

HEIESI1IS0I3- This is a broad P eak at an energy 
transfer of a few meV, where simple crystals have only 
sound waves. Glasses seem to have a sizeable amount 
of excess vibrations at this boson peak. At present, it 
is not clear which driving force brings these vibrations 
down into the low-frequency region. 

In this paper, we present for the first time experimen- 
tal evidence for a subtle but important driving force ne- 
glected so far, the vibrational entropy. 

The glass freezes at a relatively high temperature (the 
glass temperature T g ). At this temperature, the thermal 
energy is an order of magnitude higher than the vibra- 
tional energy quantum at the boson peak. Therefore any 
theoretical consideration of the origin of the boson peak 
should take the vibrational entropy into account, because 
it tends to bring the eigenmodes down to low frequen- 
cies. A lower frequency allows to distribute the thermal 
energy between more energy levels and thus increases the 
entropy 0, Ell- 
in order to see this influence in experimental data, one 
plots the eigenvalue density p(X) against the eigenvalue A 
on a logarithmic scale. Since the neutron measurement 
supplies the frequency as energy transfer E = ftw, we 
define the eigenvalue A = E 2 and measure it in meV 2 . 
The eigenvalue density p(X) — g(E)/2E, where g{E) is 
the conventional vibrational density of states. The boson 
peak is a broad peak in g{E)/E 2 . Therefore it is more 
like a shoulder in p(X). 

All measurements presented were done on the time- 
of-flight spectrometer IN6 at the High Flux Reactor of 
the Institut Laue-Langevin in Grenoble, France. The 
wavelength of the incoming neutrons was 4.1 A. Some 



data have already been published, but not in the scaling 
presented here. The samples had a scattering probability 
of about 10 % to provide a reasonable balance between 
signal strength and multiple scattering contamination. 

The data were evaluated in a new scheme for the de- 
termination of the vibrational density of states from co- 
herent and incoherent inelastic neutron scattering data, 
developed by two of us at the Institut Laue-Langevin at 
Grenoble [ll|. 

Fig. 1 shows the classical case of silica, the glass spec- 
trum of a measurement at 155 K ^2 an d the crystal spec- 
trum from a lattice dynamical calculation . In the re- 
gion around 10 meV, the crystal shows pronounced van- 
Hove singularities. The disorder should broaden these 
low-lying van-Hove singularities into a quasi-constant 
eigenvalue density on the basis of Wigner's famous so- 
lution 14] of the random-matrix problem. Instead, we 
observe a slow decrease towards higher eigenvalues, lin- 
ear on this logarithmic scale over nearly a decade. Our 
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FIG. 1: Eigenvalue density in crystalline a-quartz [l3jl and in 
vitreous silica [l2|. The line is a fit in terms of cq. @. 
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results are in qualitative agreement with recent nuclear 
inelastic scattering data [la ], which find a stronger than 
l/.E'-dependence in g(E)/E 2 above the boson peak. 

Let us see whether we can justify such a linear decrease 
on the basis of the concept of vibrational entropy. The 
tendency to form a constant eigenvalue density from the 
randomness of the dynamical matrix can be described in 
terms of an energetic level repulsion term [2J, Il4| . We 
assume this level repulsion term to be proportional to 
the eigenvalue density, with a proportionality factor B. 
The free energy Fx of an eigenmode at the eigenvalue A 
is then 

F X = Bp(\) + ^- In A. (1) 

The second term is the vibrational entropy 0, fiol ] in the 
classical limit E « k B T . 

From eq. we obtain the eigenvalue density 

which is the linear relation to In A we are looking for. 

The slope —kBT/2B should be proportional to the 
temperature, at least as long as the eigenvalue density is 



able to adapt itself to a new thermodynamic equilibrium. 
This can be used to check the validity of the scheme in 
the undercooled liquid. 

Fig. 2 (a) demonstrates the strong temperature change 
of the eigenvalue density in the heavily studied example 
of 1,4-polybutadiene [la]- The temperature dependence 
of the fitted slope k B T/2B in Fig. 2 (b) follows the ex- 
pectation of eq. (J2J) over a wide temperature range, from 
high temperatures in the liquid down until the sample 
freezes at the glass temperature. 

Polybutadiene is a clearcut example for the validity of 
eq. over a wide temperature range, with a level repul- 
sion constant B which does not depend on temperature. 
This is not always the case. One finds a slight tempera- 
ture dependence of B in the third example, selenium, as 
shown in Fig. 3. These are already published data (17j |. 
complemented by new experiments with a thinner sam- 
ple. In this case, the proportionality of the fitted slope 
to temperature is limited to a small temperature range 
around T g , extending to slightly lower values as T g . But 
within that range, the selenium data again support the 
physical picture of a balance between vibrational entropy 
and level repulsion at frequencies above the boson peak. 

If this picture is true, the boson peak must be a 




FIG. 2: (a) Eigenvalue density in crystalline, glassy and liq- 
uid polybutadiene [l6j|. Lines are fits in terms of eq. (|SJ. 
(b) Temperature dependence of the fitted slope ksT /2B in 
polybutadiene. 



FIG. 3: (a) Eigenvalue density in crystalline, glassy and liquid 
selenium |17|l . Lines are fits in terms of eq. @. (b) Temper- 
ature dependence of the fitted slope fcsT /2B in selenium. 
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crossover phenomenon between such a balance at higher 
frequency and a quite different physical reality at low fre- 
quency. We know from many experiments that we have 
well-defined sound waves at frequencies below the boson 
peak. The question is: What is the physical mechanism 
which causes such a crossover? Is it the anharmonicity 
which reduces the number of levels of a given eigenmode 
potential as one approaches the eigenvalue zero? Or is 
it the interaction between resonant modes and sound 
waves? 

Our experimental evidence does not provide a direct 
answer to this central question. But it suggests to look 
for a reasonable interpolation scheme between low and 
high frequency. At low frequency, one has again a reliable 
theoretical picture, namely an elastic medium with well- 
defined sound waves coexisting with a small number of 
additional modes. There, the neutron spectra should be 
well described by the soft-potential expression 
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This expression for the low- frequency spectrum is de- 
rived from the soft-potential model |lS|. an extension of 
the tunneling model of the low temperature anomalies of 
glasses. The first term contains the Debye frequency u>d 
and describes the sound waves. In addition to the sound 
waves, the model postulates a continuous distribution of 
additional modes around the eigenvalue zero. The posi- 
tive eigenvalues provide vibrational resonant modes coex- 
isting with the sound waves, the second term of eq. J3J • 
The negative eigenvalues, in principle unstable modes, 
are supposed to be stabilized by the anharmonic fourth 
order term of the mode potential. They lead to double- 
well potentials, giving rise to tunneling states at low tem- 
peratures and to classical relaxation at higher tempera- 
tures. The third term of eq. (J3J) is the soft-potential 
expectation for this classical relaxation spectrum. f v u, 
and f re i are given in terms of the parameters of the soft- 
potential model |is| 
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Here P s is the density of additional modes around the 
eigenvalue zero, M is the average atomic mass, p is the 
mass density and W is the crossover energy between vi- 
brational and tunneling states at low temperatures. 

The soft-potential model has been checked against the 
low-tempcrature glass anomalies in the specific heat, the 
thermal conductivity and in the mechanical loss 0, 0] . 
The model predictions were found to be essentially cor- 
rect, with an important exception: as soon as the barrier 
height of the double-well potentials begins to be a size- 
able fraction of the thermal energy at the glass transition, 



the measured classical relaxation gets much weaker than 
the soft-potential prediction . As we will see, the 

same effect appears in the neutron spectra at elevated 
temperature. 

In order to check the neutron data, one needs a suitable 
interpolation scheme between eq. Q below the boson 
peak and eq. (J2J above. As it turns out, it is convenient 
to use 
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where PhighW is given by eq. ©. 

Eq. © is an elegant way of interpolating from one 
eigenvalue density to the other, without introducing an 
additional crossover parameter. One can characterize the 
whole spectrum with five reasonable parameters, each of 
which has a well-defined physical meaning either in the 
low-frequency or in the high-frequency range. 

The above-mentioned failure of the soft-potential 
model is most clearly seen in the neighborhood of the 
glass temperature. Fig. 4 shows the neutron eigenvalue 
spectrum of B2O3 at its glass temperature of 550 K [2lj . 
together with the sound wave expectation on the basis 
of Brillouin scattering measurements at the same tem- 
perature [23I . There is a sizeable relaxational signal at 
the low-frequency end, which allows to fit / re ; with an 
accuracy of 25 %. The vibrational contribution f v u, is 
fitted with the same accuracy. If one compares the ratio 
f rel/ fvib with the soft-potential expectation of eqs. (@J 
and 10, it is a factor of five too low (the comparison re- 
quires the crossover energy W between vibrational and 
tunneling states, which was taken from ref. 19]). The 
same is true for polystyrene at its glass temperature 360 
K [23| in Fig. 4 and for the case of 1,4-polybutadiene at 
its glass temperature of 180 K in Fig. 2 (a). 
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FIG. 4: Eigenvalue density in B2O3 12 II and in polystyrene 
(2^| at the glass transition. Lines are fits in terms of eq. 
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While these results cannot be understood in terms of 
the soft-potential model, they support independently our 
conclusion on the importance of the vibrational entropy. 
f v ib describes the low-frequency additional vibrations at 
low energy transfer E. These are eigenmodes which still 
have a positive eigenvalue. Therefore they profit in full 
from the vibrational entropy, having an entropy contri- 
bution to the free energy of ksT\n.(E /fc^T). This is 
different for f re i, which stems from eigenmodes with a 
negative eigenvalue. These modes correspond to relax- 
ational jumps between neighbouring energy minima, with 
an entropy contribution to the free energy of the order 
— fcsTln2. Thus one expects a Boltzmann factor ratio 
of the order of kBT g /2E between positive and negative 
eigenvalues close to zero, if both freeze under otherwise 
identical conditions at the glass temperature. Since the 
thermal energy at the glass transition is an order of mag- 
nitude larger than the energy transfer at the boson peak, 
this is the observed factor of five. 

To summarize, we find threefold experimental evidence 
for the influence of vibrational entropy on the neutron 
spectra of glass formers. The first is a decrease of the 
eigenvalue density above the boson peak, linear in a plot 



against the logarithm of the eigenvalue (this has already 
been seen in nuclear inelastic scattering [Lsl ]. though fit- 
ted differently). The second is the temperature depen- 
dence of the corresponding slope above the glass temper- 
ature, proportional to temperature in 1,4-polybutadiene 
and selenium. The third is the weakness of the relax- 
ational part of the spectrum close to the glass transition, 
which is a factor five weaker than the soft-potential ex- 
pectation of an equal footing of positive and negative 
eigenvalues close to the eigenvalue zero. 

These findings do not resolve the boson peak problem, 
but they set the stage for a proper theoretical treatment. 
They show the existence of a theoretically well-defined re- 
gion at frequencies above the boson peak, where crystals 
have their first van-Hove singularity. At the glass tran- 
sition, this region is dominated by the balance between 
level repulsion and vibrational entropy; sound waves are 
heavily overdamped there |24|. On the other hand, the 
sound wave picture is the correct starting point in a glass 
below the boson peak frequency, where the glass is an 
elastic medium with a small number of additional exci- 
tations. The boson peak must be the crossover between 
these two regions. 
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